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ABSTRACT 

We present the results of a systematic search for transiting planets in a ~5 Myr 
open cluster, NGC 2362. We observed ~1200 candidate cluster members, of which 
~475 are believed to be genuine cluster members, for a total of ~100 hours. We identify 
15 light curves with reductions in flux that pass all our detection criteria, and 6 of 
the candidates have occultation depths compatible with a planetary companion. The 
variability in these six light curves would require very large planets to reproduce the 
observed transit depth. If we assume that none of our candidates are in fact planets 
then we can place upper limits on the fraction of stars with hot Jupiters (HJs) in 
NGC 2362. We obtain 99% confidence upper limits of 0.22 and 0.70 on the fraction of 
stars with HJs (f p ) for 1-3 and 3-10 day orbits, respectively, assuming all HJs have a 
planetary radius of 1.5i?j up . These upper limits represent observational constraints on 
the number of stars with Hjs at an age <10 Myr, when the vast majority of stars are 
thought to have lost their protoplanetary discs. Finally, we extend our results to the 
entire Monitor Project, a survey searching young, open clusters for planetary transits, 
and find that the survey as currently designed should be capable of placing upper 
limits on f p near the observed values of f p in the solar neighbourhood. 
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1 INTRODUCTION 

1.1 Transiting Planets Around a Young Star 

The detection of a transiting planet provides a wealth of 
information which cannot be matched by any other plan- 
etary detection method at the moment. There have been 
over 20 transiting planets detected thus far, and the list 
is growing rapidlyj. These systems provide the only means 
for directly measuring a planet's radius. They also pro- 
vide a means for measuring the inclination of the orbital 
plane, which in turn removes the sin i ambiguity of a planet 
only detected via the radial velocity method. If the host 
star is bright enough, transiting planets also provide a 
means for measuring their atmospheric composition when 
the stellar spectrum, in eclipse (planet behind the star), 
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is subtracted from the combined- li ght out of eclipse spec- 
trum of the star plu s the planet dRichardson et ahl 120071 ; 
iGrillmair et al. 1 120071 1 . Additionally, very precise measure- 
ments of the time of transit can be used to search for 
and characterise the orbits of other, sometimes very low- 
mass, planets in the extrasolar system, because these other 
planets would slightly perturb the transiti ng planet's orbit 
i|Agol et al.ll2005l : lHolman fc Murravll2005D . 

Despite the recent success of many transit surveys, there 
has yet to be an observation of a transiting planet orbiting 
a PMS hydrogen burning star. In fact, there is a relative 
paucity of detected planets orbiting any stars at an age less 
than ~100 Myr. Aside from the re cent claim of the detec tion 
of a planet orbiting TW Hydrae (|Setiawan et al.ll200St ). we 
are unaware of any detections of a planet orbiting a <100 
Myr hydrogen burning star. 

Planets are believed to form in protoplanetary discs 
where formation and accretion, as well as migratio n, halt 
following the dispersal of the disc. lHaisch et alTi|200ll) found 
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that roughly 50% of stars lose their discs by ~3 Myr, while 
nearly all discs are dissipated by ~10 Myr. This places a sig- 
nificant constraint o n the formation time of a gaseous giant 
planet at <10 Myr feodenheimer fc Linll2002l '). The discov- 
ery of a planet around a very young star, specifically ~5 
Myr, would provide important constraints for planetary for- 
mation mechanisms, migration time scales and dynamical 
evolution, a nd their relation to disc l ifetimes and clearing 
time scales (|Bodenheimer fc Linl 12002'). A transiting planet 
at ~5 Myr would lead to constraints on the planet mass and 
radius at an early age, i.e. very near the initial conditions for 
giant planet evolution, which are essentially unconstrained 
at the moment. 



1.2 The Cluster 

NGC 2362 is a well studied open cluster wh ose age (~5 Myr; 
iMoitinho et al.ll200ll . iDelgado et alj|2006to . which coincides 
with the tail of the distribution of circumstellar disc life- 
times, a nd relatively moderate distance (1480 pc, (m — M)q 
= 10.85; IMoitinho et al.ll200lh make it an ideal test bed for 
the study of stars and their environments during the early 
pre-main sequence (PMS). For example, the detection of de- 
tached eclipsing binary (EB) stars in the cluster will provide 
measurements of each star's mass and radius without the 
use of a model. These measurements can, in turn, be used 
to constrain PMS evolution models. 

Many recent surveys have been conducted to charac- 
terise the properties of NGC 2362. The most relevant work 
includes the determin ation of fundamental cluster properties 
l|Moitinho et al.ll200ll ). the determin ation of the circum stel- 
lar disc fraction from near-IR excess (|Haisch et al.ll200ll ), an 
Ha e mission surv ey to study T Tauri stars and disc accre- 
tion (|Dahmll2005l l. and a study of primordial circumstellar 
discs usin g infrared excesses measured by the Spitzer Space 
Telescope (|Dahm fc Hillenbrand|[2007h . 

1.3 The Survey 

We have completed a high cadence photometric monitoring 
survey of NGC 2362, with observations made using the Mo- 
saic II imager on the 4m Blanco telescope at CTIO. There 
are three primary scientific goals of this survey: one, to dis- 
cover low-mass EB systems which will allow us to simulta- 
neously measure the mass and radius for each member of 
the system, two, to search for transiting planets orbiting 
PMS stars, and three, to characteris e the rotation pe riods 
for low-mass members of the cluster (jlrwin et al.l 2008). 

Our observations of NGC 2362 were designed to be sen- 
sitive enough to detect a large planet (R p > 1- Rj up ) in a 'very 
hot- J upiter' orbit, i.e. orbital period < 3 days (|Aigrain et al.l 
I2007T ). Specifically, given the average observational cadence 
of ~6 min during 18 nights spread over a year with opti- 
mal signal-to-noise we ought to be able to photometrically 
detect planets with R p > lRj up aroun d a 0.7 M^ primary . 
For lower mass, and hence fainter, stars lAigrain et all (120071 ) 
predicted that we should be able to detect planets with 
-Rp > 2i?j up orbiting a 0.2 Mq primary. This corresponds 
roughl y with the spectrosc opic limits from 8 m class tele- 
scopes. lAigrain et al.l l|2007h shows that it would be possible 
to detect the radial velocity (RV) signal of a lMj up planet 



orbiting a 0.2 Mq primary in a 1 day orbit down to I ~ 
18 using the UVES spectro graph on the Very L arge Tele- 
scope. Initial simulations bv lAigrain et alj (|2007l ) predicted 
that there would be 4.7 EBs and zero transiting planets in 
our data. 

These observations are part of a larger photometric 
survey of 9 young (< 200 Myr) open clusters covering a 
wide range of ages a nd metallicities (the M onitor projec10; 
iHodgkin etal]|2006l and lAigrain et al.ll2007l ). 



1.4 Methods 

We developed an automated method to search for occul- 
tations in the light curves generated from our NGC 2362 
observations^ In particular, we attempt to remove the sig- 
nificant spot-induced variability displayed by many of our 
young and late-type targets. As transits and eclipses are es- 
sentially the same from a detection standpoint, and cannot 
always be distinguished from the light curve alone, we de- 
tect both types of events among our candidates. We defer 
the discussion of the EB candidates to a later paper. Here 
we focus on our transit candidates, for which spectroscopic 
follow-up observations, to confirm cluster membership and 
ascertain the source of the occultations by measuring the 
mass of the occulting body, are underway. 

We then perform extensive Monte Carlo simulations to 
evaluate our sensitivity to planets of different radii and pe- 
riods around stars of different masses in the cluster, and 
use these simulations to place upper limits on the incidence 
of hot and very hot Jupiters (HJs and VHJs) at 5 Myr. 
These upper limits represent constraints on HJ incidence at 
an age of less than 10 Myr, when any HJs should have re- 
cently finished migrating toward their host star. They there- 
fore constitute an important measurement for constraining 
planetary formation and migration time-scales. 



1.5 Organisation of Paper 

The remainder of this paper is structured as follows: the 
data reduction process is discussed in Section 2, and our 
method of cluster membership selection is presented in Sec- 
tion 3. The technique used to search for occulting systems 
and the resulting candidates are described in Section 4. The 
Monte Carlo simulations and the derived upper limits on the 
fraction of stars in NGC 2362 with short period planets, are 
discussed in Section 5. Here we also compare our sensitivity 
and results to other transit surveys that have targets stellar 
clusters. Finally, we present our conclusions in Section 6. 
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3 From this point forward, when we refer to eclipses we specifi- 
cally mean the eclipse of one star by another star or brown dwarf, 
while transits only refer to the case of a planet transiting a star. 
We use the term occultation in a more general sense to include 
both of these phenomena, but to also include all situations where 
the flux from a star has been reduced because some other body 
has passed in front of the stellar disc. 



Transits in NGC 2362 3 



Night 


Frames / night 


Start 


Finish 


1 


50 


406.03573 


406.29499 


2 


47 


407.03439 


407.28675 


3 


42 


408.03985 


408.30530 


11 


35 


416.04713 


416.24360 


12 


41 


417.02969 


417.24706 


13 


42 


418.04069 


418.24918 


323 


24 


728.21098 


728.35062 


326 


30 


731.21755 


731.36089 


328 


20 


733.27306 


733.36723 


330 


30 


735.20998 


735.36427 


331 


32 


736.20901 


736.36480 


332 


29 


737.21096 


737.36434 


333 


31 


738.21689 


738.36825 


336 


31 


741.21493 


741.36774 


359 


62 


764.04498 


764.33686 


360 


19 


765.03671 


765.33456 


361 


61 


766.03490 


766.32738 


362 


72 


767.03074 


767.33916 



Table 1. Summary of photometric observations of NGC 2362. 
Night is the night of observations relative to the first night of 
observations. Start and Finish are the beginning and end of the 
nightly observations, respectively, given in HJD - 2453000.5. 



2 OBSERVATIONS AND DATA REDUCTION 

Photometric monitoring data were obtained using the CTIO 
4m Blanco telescope, with the Mosaic-II imager, during 18 
nights from February 2005 to January 2006. This instrument 
provides a field of view of ~36'x 36' (0.37 deg 2 ), using a 
mosaic of eight 2k x 4k pixel CCDs, at a scale of ~0'.'27 per 
pixel. A summary of our observations, including the number 
of frames per night as well as the start and finish time for 
each night, is given in Table Q] 

Our observations consist of a series of 75 s i-band expo- 
sures with an average cadence of ~6 min. We also obtained 
a few longer images in the V-band (3 x 600s) in photometric 
conditions for the production of a colour magnitude diagram 
(CMD). Our observations are sufficient to give 1% or better 
photometric precision per data point from saturation at i ~ 
15 down to i ~ 19, as seen in Figure [T] This range corre- 
sponds to G through mid-M spectral types at the age and 
distance of NGC 2362. 

For a full descript ion of our data reduction procedure 
see llrwin et al.l l|2007aT ). Briefly, after we corrected for cross- 
talk between the detector readouts, we followed the standard 
CCD reduction scheme of bias correction, flat-fielding, and 
astrometr ic calibration before p hotometric calibration as de- 
scribed in llrwin fc Lewis! (|200lh . Following this we generated 
a master catalogue for the i-band filter by stacking 20 of the 
frames taken in the best conditions (seeing, sky brightness 
and transparency) and running the source detection software 
on the stacked image. The resulting source positions were 
used to perform aperture photometry on all of the images, 
with the final result a time-series of differential photometry. 
We achieved a per data point photometric precision of ~2-4 
mmag for the brightest objects, with RMS scatter < 1% for 
i < 19 (see Figure [TJ. 

Our source detection software flags any objects detected 
as having overlapping isophotes as likely blends. This in- 
formation is used, in conjunction with a morphological im- 




i magnitude 



Figure 1. Plot of the rms scatter per data point for the entire 
series of observations as a function of magnitude for the i-band 
observations of NGC 2362. We include only unblended objects 
with stellar morphological classifications. The diagonal dashed 
line shows the expected rms from Poisson noise in the object, the 
diagonal dash-dot line shows the expected rms from sky noise in 
the photometric aperture, and the dotted line shows an additional 
1.5 mmag contribution added in quadrature to account for sys- 
tematic effects. The solid line shows the total predicted rms from 
these effects. This plot shows the rms of our light curves before 
they have been corrected for seeing. 



age classification flag also generated by the pipeline software 
i Irwin fc Lewis! l200lh . to allow us to identify non-stellar or 
blended objects in the time-series photometry. 

The CCD magnitudes were converted to the stan- 
dard Johnson- Cousins system using regular observations of 
lLandoltl (|l992t ) equatorial star fields in the usual way. 

Light curves were extracted from ~85 000 objects, 56 
000 of which had stellar morphological classific ation, using 
our st andard aperture photometry techniques (jlrwin et alJ 
12007 j ). We fit a 2D quadratic polynomial to the residuals 
in each frame (measured for each object as the difference in 
magnitude between the current frame and the median taken 
over all the frames) as a function of position, for each of 
the eight CCDs separately. We then removed this function 
to account for variations in transparency and differential 
atmospheric extinction across each frame. For a single CCD, 
the spatially varying part of the correction remains small, 
typically ~0.02 mag peak-to-peak. 

As a last step there is a small correction applied to all 
the light curves for seeing-correlated effects. This was done 
by looking for seeing-correlated shifts in the light curve from 
its median magnitude. A simple quadratic polynomial was 
fit to the shift as a function of the full width half max of the 
stellar images on the corresponding frame. This fit was then 
subtracted from the light curve. Typically, this fit would 
reduce the rms of the light curve by < 0.01 mag, however, for 
the cases that showed the strongest correlations with seeing 
the reductions in rms would be fairly significant, ~0.1 mag. 

For the production of deep CMDs, we stacked 20 i-band 
images taken in good seeing and photometric conditions. 
The long exposure T/-band frames were stacked before run- 
ning source detection and the accompanying astrometric and 
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photometric calibration. The l/-band sources were matched 
against sources in the i-band master catalogue, which then 
enabled us to produce a deep CMD. The limiting magni- 
tudes, measured as the approximate magnitude at which 
our c atalogues are 50% complete were V ~ 24.4 and i ~ 
23.6 i|lrwin et al.ll2007al ). 



3 SELECTION OF LOW-MASS CANDIDATE 
MEMBERS 

Before we could search for transits we had to iden- 
tify low-mass cluster members. List s of candidate mem- 
bers are available in the l i terature |Moitinho et al.l l200ll . 
iDahm fc Hillenbrand! 120071 . iDeleado et alJ l2006h . however, 
in order to match the field of view of our survey, which is 
wider than previous surveys of NGC 2362, we elect to use a 
V versus V — I CMD for candidate membership selection. 



3.1 The V versus V-I CMD 

The V, V — I CMD used for candidate membership selection 
is shown in the upper panel of Figure [2] The V and i mea- 
surements were converted to the standard Johnson- Cousins 
photometric system using colour equations from our obser- 
vations of photo metric standard stars (see Eqns. 1-3 from 
llrwin et alj|2008l ). 

The cluster sequence is difficult t o identify by eye , espe - 
ciall y when compa red to the CMDs of lMoitinho et all (|200ll ) 
and lPahml (2005). Our survey covered a much larger por- 
tion of the sky (0.3 7 deg. 2 compared t o ~0.05 deg . 2 for 
iMoitinho et~al]|200ll and ~0.03 deg. 2 for |Pahm|[2(5o5l ') and 
therefore suffers from greater field contamination. Thus, we 
used a second CMD including only objects within 7 arcmin 
of t CMa, defined as the centre of the cluster, to identify 
the cluster sequence, shown in the lower panel of Figure [2] 
The 7 arcmin radius was chosen based on a 2a cut from 
a Gaussian fit to the radial profile of rotating stars in the 
field of view. Short period rotators, which are indicative of 
youth and therefore also c luster membership , were selected 
from the identifications by llrwin et alJ (|200ST ) . The presence 
of rotators at large (> 15 arcmin) distances from the cen- 
tre of the cluster justifies the use of the entire field of view 
in our candidate selection process, despite the ensuing high 
contamination from field objects. We elected to do this be- 
cause we did not want to miss any transiting planets, which 
are rare in both the cluster and the field. 

We t hen followed the c andidate selection method de- 
scribed in llrwin et a l. (2008). Briefly, we manually defined 
an empirical cluster sequence that follows the sequence visi- 
ble in Figure[2jD. All objects falling between two cuts defined 
by shifting the sequence right and left (perpendicular to the 
sequence itself) were then selected as candidate members. 
This led to the selection of 1813 candidate members over 
the full range from V=15.7 to 26. 

To determine model masses and radii for our can- 
didates, the I-band absolute magnitudes of the NextGen 
model are used, because these are less susceptible to a miss- 
ing source of opacity, which creates a discrepancy between 
the models and observations in the V — I colour for T a ff 
< 3700 K (corresponding in this case to V — I > 2.5 
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Figure 2. (a) V versus V-I CMD of NGC 2362 from stacked 
images for all objects with stellar morphological classification. 
The cluster boundaries, which define candidate membership, are 
shown as dashed lines (all objects between the lines were selected). 
The cluster sequence is clearly contaminated by objects in the 
galactic field, (b) V versus V—I CMD of objects within ~7 arcmin 
of r CMa, here defined as the centre of the cluster. The 7 arcmin 
cut was determined by the spatial distribution of rotating stars, 
as described in the text. The cluster sequence is clearly visible 
and our empirical sequence is shown as the solid line. The mass 
scale in both plots is from the 5-Myr NextGen model isochrone 
llBaraffe et aljl 998l, using our empirical isochrone to convert the 

V magnitudes to I magnitudes, and subsequently obtaining the 
masses from the I magnitudes due to known problems with the 

V magnitudes of the models (see Section 3.1). 
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IBaraffe et al.lll998h . Therefore the /-band absolute magni- 
tudes give the most robust estimates of mass and radius. 
Our adopted mass-radius -magnitude relation comes from 
the procedure described in lAigrain et al. I d2007t) . which com- 
bines the NextGen isochro nes of Baraffe et all rtl9 98). with 
the DUSTY isochro nes of IChabr er et all (|200Ch . and the 
COND isochrones of IBaraffe et all (|2003l ~This adopted re- 
lation covers a mass range from 0.5Mj up to 1.4Mq. For the 
few, ~30, objects that are brighter than the limits of the 
NextGen models (i.e. Mi < 3.9, o r M» > 1.4M ) we use 
the isochrones of ISiess et al.l ((2000) to determine the model 
masses and radii of those particular candidate cluster mem- 
bers. 



3.2 Contamination 



I Irwin et al.l (|2008T ) estimate the level of contamination for 
the sample of candidate cluster members to be ~65%. They 
caution that their estimate is s omewhat un c ertain due to the 
need to use Galactic models. Ilrwin et al. also note 

that spectroscopic follow-up will be needed to make a more 
accurate contamination estimate. We note that when the 
range of magnitudes is restricted to those that we search 
for occultations (see Section 4) the contamination level is 
slightly reduced to ~60%. 



4 OCCULTATION DETECTION 

While we are nominally searching for transits, our search 
procedure identifies any light curve with occulting events. 
Therefore, throughout this section we will discuss our search 
for occultations, which will inevitably yield a list of transit 
candidates. After achieving the necessary signal-to-noise, re- 
moving systematic trends, and obtaining sufficient coverage 
and sampling in the data, perhaps the most significant ob- 
stacle in any systematic search for occulting systems is the 
intrinsic variability present in many stars' light curves. Peri- 
odic variability, typically due to the rotation of spots on the 
surface of the star, is a particularly severe contaminant be- 
cause it leads to regular reductions in the observed flux from 
the star, which is precisely the behaviour (i.e. occultations) 
we are trying to identify. Non-periodic variability is also sig- 
nificant in young stars: for instance, a star could change 
brightness following occ ultations by or intera ctions with cir- 
cumstellar material (see lBouvier et aL I l2007l for a discussion 
of AA Tau, a young star whose photo metric variability orig- 
inates from interactions with its disci iDahm fc Hillenbrand] 
l|2007l ) found an upper limit of ~7% for t he fraction of star s 
with optically thick discs in NGC 2362. Ilrwin et all (|2008h 
note the difficulties in determining the fraction of cluster 
members that rotate, because the rotation sample has a 
lower contamination level than the remainder of the general 
candidate cluster members. Applying the same contamina- 
tion estimates to both populations they estimate a conserva- 
tive lower limit of ~14% for th e fract ion of cluster members 
which are rotating. Ilrwin et al. I (|2008l ) estimates that the ac- 
tual fraction of cluster members which rotate is ~40% based 
on the high correlation between rotation and cluster mem- 
bership. Clearly, there are more rotators than stars with 
optically thick discs, which is fortunate because rotational 
variability is easier to filter than non-periodic variability. 



The effects of rotation lead to smooth variations in bright- 
ness. When this is the dominant source of variability it can 
be subtracted from the signal without introducing signifi- 
cant additional features into the light curve. 

Before we began our search we removed from our sam- 
ple a number of lig ht curves which wer e flagged by the data 
reduction process (|lrwin et all l2007al ). Specifically, we re- 
moved any light curves that were flagged as saturated and 
any systems where more than 10% of the data points be- 
longed to low confidence regions in the standar d size aper- 
ture a s flagged by our reduction procedure jlrwin et all 
2007a). We found a few stars which displayed behaviour 
consistent with saturation were not flagged in the original 
procedure. Thus, we elected to visually examine the remain- 
ing light curves by eye and flag those that were saturated. 
This led to the removal of an additional 75 objects. Objects 
that were flagged as blended were not excluded, how ever, 
because the reduction procedure of Ilrwin et all (|2007al ) did 
a sufficient job in correcting the effects of blending such that 
these objects could be searched for occultations along with 
the rest of the sample. 

We also limited our search with a magnitude cut such 
that we only examine objects brighter than / = 19. We place 
these cuts based on the limitations of spectroscopic follow- 
up: the cause of occultations must be confirmed with RV 
measurements (for a full des cription of the spec troscopic 
limitations of this study see lAigrain et all 120071 ) . Briefly, 
lAigrain et ail (|2007T ) found that with existing multi-object 
spectrographs on 8 m class telescopes it would be possi- 
ble to reach RV precision of ~ 2 km/s down to / ~ 18. 
Our cut of I = 19 is thus conservative, yet we point out 
that ISahu et all (|2006h were able to measure RV variations 
to ~1 km/s for an object with / = 18.75 using the UVES 
echelle spectrog r aph at the 8-m Very Large Telescope, while 
IWeldrake et all (|2008al ) were able to detect a K = 114 m/s 
planet in a 4-day orbit around a K star of V = 17.4. We note 
that any candidates fainter than ~18 in I will be extremely 
difficult to follow-up in order to confirm the origins of the 
occulting behaviour found in the light curve. 

Following the removal of these objects there remained 
1180 candidate cluster members to be searched for occulta- 
tions. 



4.1 Variability Filtering 

The difficulty with attempting to remove the intrinsic stellar 
variability from a light curve is that many filtering methods 
will remove or affect the occultation signal as well as the pho- 
tometric signal coming from just the star alone. Therefore, 
the challenge lies in developing a method that successfully 
disentangles the occultation signal from that of the host star. 

Because we find rotation to be the dominant source of 
variability in the candidate members of NGC 2362, we fo- 
cused on removing this source of variability. After ka clip- 
ping outlying data points in the light curves, we selected 
rotating stars in our sample by performing a least-squares 
sine-fit to the time series m(t) (in magnitudes) of every can- 
didate cluster member using 



m(t) = rridc + a sin(2nt/P + (f>), 



(1) 



where rridc is the mean light curve level, a the amplitude, 
<j> the phase and P the period of rotation. For the fits rridc, 
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5-2085 P = 2.941 a, = 0.018 



<j 2 = 0.027 




-M 1 



1000.0 




Figure 3. Example rotation fit to star 5-2085. Shown is the light 
curve and best fit folded over the rotation period (top) and the 
period-folded sine-subtracted light curve (bottom). The three dif- 
ferent observing windows, as described in the text, are shown with 
0.15 mag offsets in the top panel. It is clear to see that both the 
amplitude and phase of the signal from this star change with time. 
P is the best fit period in days while on is the best fit amplitude 
in magnitudes for each of the observing windows. 



a, and <fi are & ee parameters at each value of P over a grid 
of equal logarithmically spaced steps in period from 0.1-181 
d (corresponding to half the time between our first and last 
observations). Our fits adopted a single period, but we al- 
lowed the phase and amplitude to change following any gaps 
in our observations of more than three weeks. We allowed 
these changes because the size and location of star spots can 
evolve very rapidly over these time scales in young stars. We 
fix the period, however, because we would not expect a sig- 
nificant change in the angular momentum of the star in the 
course of a single year. The output of this procedure is a 
'least-squares periodogram,' and the best-fitting period is 
the one with the lowest reduced \ 2 ■ An example of this fit- 
ting procedure is shown in Figure [3] where the data have 
been folded on the best fit period. From the figure it is easy 
to see the change in amplitude and phase following long gaps 
in our observations. 

We measured the reduced \ 2 m our light curves before 
(x^.flat) and after (xtm) we subtracted the best sine-fit and 
selected rotators based on the change in reduced % 2 '- 

Ax'/xLat > 0.7, (2) 

where x^.flat is the reduced % 2 °f the original light curve 
with respect to a constant model, and A\t is the change in 
reduced x 2 following the fit. We also required that xS.at < 
60 for an object to be classified as a rotator. We acknowl- 
edge that a reduced \ 2 °f 60 is quite large, however, we 
find that these systems (xt at ~ 60) display periodic vari- 
ability. The large value of \ 2 ls the result of extremely high 
signal-to-noise data and slight departures from the mod- 
elled behaviour of Eqn. [T] The upper limit for xt.at was 
selected because we found that Eqn. [T] was a poor model for 
all objects above this threshold. Figure [3] shows Ax 2 plotted 
against xt at > an d highlights the systems which were selected 
as rotators. 

The empirically determined cutoff in Eqn. [2] is more 
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Figure 4. as a function x 2 , at ^ or a 'l candidate cluster mem- 

bers in our sample. The stars highlight the systems that were 
selected as rotators. 



stringe nt than the initial cut used in the M oni tor rotation 
pape rs (jlrwin et al ]|2006l .l Irwin et al . 2007b, and llrwin et al.l 
l200St) . Those papers employ visual inspection of each light 
curve to remove any objects without clear periodic variabil- 
ity from their sample. Therefore, their samples are more 
complete. We wish to avoid visual inspection, however, in 
order to remove hu man interact i on fro m our selection pro- 
cedure. As noted bv lBurke et alj ([2006) , a clear set of detec- 
tion criteria that do not rely on human input are extremely 
important for establishing the actual sensitivity of a survey 
to planetary transits using Monte Carlo simulations. There- 
fore we chose not to add a step with visual inspection. We 
do find that there is ge nerally good agreem ent between our 
sample and the one in Irwin et al 1 20081 ). except for the 
slowest rotators. Ilrwin et al.l 1 200S ) use data from only a 
single observing season, so they are not sensitive to periods 
> 24 days. We also note that our classification scheme leads 
to some objects, whose variability is not the result of star 
spots, to be misclassified as rotators, as can be seen by the 
slight build up around 0.5 and 1 days in the period distribu- 
tion of Figure [S] This build up is caused by low amplitude 
night edge effects. For the purposes of this study this less 
than perfect classification scheme is acceptable because the 
subtraction of the sine-fits was designed to remove periodic 
variability of a non-eclipsing nature. 

Our goal was to identify any objects which clearly dis- 
played periodic variation of a relatively long temporal signa- 
ture while excluding those with high-frequency events where 
the reduced x 2 would be significantly improved by a sine-fit 
(i.e. objects with multiple eclipses). We wanted to test that 
our classification scheme did not identify transits as rota- 
tion. To do this we took a sample of our flattest light curves, 
selected for their low \ 2 (x^flat < 1-5) and low dispersion 
in median flux for each observing se ason, and inserted tran - 
sits according to the formalism of iMandel fc Ago] d2002h 
using the limb-darkening coefficients of lClaretl I 2000l) . The 



planet radius and orbital period were chosen randomly over 
the same grid as that used in the simulations in Section 5. 
The orbital inclination and epoch were chosen such that we 
would have observed at least a portion of one transit. We 
then ran the rotation test described above and we found 
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Figure 5. Rotation period versus mass for the candidate cluster 
members classified as rotators. The blue crosses show those ob- 
jects classified as rotators in th is paper. Th e large black points 
show rotators as identified by llrwin et alj fcOOSt) . |l rwin et al,l 
(2008) only use a portion of the data, so they are not sensitive 
to periods greater than ~24 days. We note that not every cross 
represents a rotator, especially the excess of objects at ~0.5 and 
~1 d (due to night-edge effects). Other discrepancies between the 
two samples are likely the result of different windows of observa- 
tion, such that in some sta rs spot activity wa s likely low during 
the subset of the data that llrwin et al. | J2008h examined. 



these objects to be misclassified as rotators in < 0.1% of 
our test cases. This extremely low rate of false classification 
allowed us to proceed with confidence that we were not re- 
moving occultation signals from our light curves. We were 
not worried about misclassifying EBs because we found that 
the removal of a sine curve from an EB light curve left large 
residuals which were detected by our occultation search al- 
gorithm. 

Following this procedure we detected rotation in 268 
stars, or roughly 23% of the candidate cluster members in 
our sample. In the objects where rotation was detected we 
subtracted the best sine-fit from the light curve and searched 
for occultations in the same manner used for non-rotators 
des cribed below. 

lAigrain fc Irwinl (|2004T l propose a number of alternative 
methods for filtering intrinsic stellar variability without re- 
moving the signal from an occultation. We attempted to use 
their least-squares filtering method and non-linear filtering 
method but we fou nd these removed the s ignal from even 
the deepest eclipses. lAigrain fc Irwinl (|2004l ) designed these 
methods for space based occultation surveys, and caution 
that their procedures are likely only valid for observations 
which are continuous on time scales 3> a typical occultation. 
Unfortunately this is rarely, if ever, the case with ground 
based surveys and we confirm their initial forewarnings. 

Lastly, we also attempt ed to model th e variations us- 
ing the formal spot model of iDorre n (1987). We found this 
method to be far too computationally intensive for an auto- 
mated procedure, while the least-squares sine- fit serves as a 
good proxy to the full spot model. 



4.2 Noise Properties 

Before we began our search for occultations we examined 
our light curves for the presence of correlated, or red, noise. 
Ground b ased surveys have been shown to suffer from red 
noise (see IPont et"al] 120061 for a very detailed discussion, 
and the typical red noise levels in several existing surveys), 
which invariably makes it more difficult to detect occulting 
objects. This correlated noise means that the uncertainty 
in data binned over n points decreases slower than in the 
uncorrelated, or white, noise case, where the uncertainty in 
binned data is oc 1/y/n. 

Our search for occultations is based on a Box Least 
Squares (BLS) fit, where the box represents a short time 
scale periodic decrease in the mean flux from the star. There- 
fore correlated noise on the same time scale as an occultation 
can lead to a large detection statistic for every light curve, 
even those which are spurious. A large detection statistic in 
every light curve necessitates an extremely large detection 
threshold, meaning that only the most significant occulta- 
tions are followed up, while shallow occultations or light 
curves with only a few o ccultation data points are not de- 
tected. [Pontj?tjil] |2006) showed that it is possible to mod- 
ify the standard white noise detection statistic to account 
for correlated noise, however, thereby eliminating many of 
the spurious candidates. 

We characterise the presence of red noise in o u r sur - 
vey according to the method outlined in lPont et al.l (2006) 
by examining the flattest light curves (low reduced \ 2 with 
respect to a constant model) in our sample. These ob- 
jects, which exhibit little variability, should be dominated 
by noise. Figure [6] shows the RMS scatter for individual 
points as a function of magnitude as well as the RMS in 
15-adjacent-point averages (which for the sampling of our 
data corresponds to roughly 2.5 hours, a typical time scale 
for transiting hot Jupiters) , compared to the expected value 
of the RMS in 15-adjacent-point averages in the presence 
of white noise. It is clear to see that the expected l/s/n de- 
crease in the noise does not apply to most of our light curves. 
Therefore, an assumption of white noise is unfounded, and 
any attempt to detect occultations must take this red noise 
into consideration. 

For our survey we find correlated noise at a level of ~1 
mmag for the brightest stars in our sample. This is on a 
similar scale to t he best ground based surveys discussed in 
IPont et al.l (|2006T l. 



4.3 Search for Occultations 

Following the removal of periodic variability we searched 
all the candidate cluster members fo r occultations usin g a 
refined version of the BLS algorithm dKovacs et al.ll2002l) as 
designed bv lAigrain fc Irwinl I (|2004h . lAigrain fc Irwinl "(2004) 



show that the best fit model reduces to finding the inverse 
variance weighted mean of the observed occultation data 
points. They test the significance of a detection using the 
familiar detection statistic S: 



(3) 
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Figure 6. RMS as a function of magnitude for a subset of candi- 
date cluster members in NGC 2362. The filled circles represent the 
RMS scatter per data point, the triangles the RMS for 15-point 
averages, and the stars represent the expected values of the 15- 
point averages assuming white noise. The triangles lie well above 
the stars indicating the presence of red noise over ~2.5 hour time 
scales. For the brightest stars this appears to be the dominant 
effect in the noise. 



where the sum includes all in-occultation data points i, di 
is the difference in flux between the ith data point and the 
mean flux of the entire light cu rve, and ai is the uncert ainty 
in the ith flux measurement. lAigrain fc Irwinl (|2004 ) also 
show that S 2 is equal to the difference in \ 2 between a flat 
model and the best fit occultation model of the data. De- 
tections with large S are typically considered the be s t can - 
didates for spectroscopic follow-up. lAigrain fc Irwinl (|2004l ) 
only considers the case of white noise, however, which can 
lead to many spurious candidates. We discuss the signifi- 
cance of a detection in the presence of red noise below. 

We restricted our search to occultations that happen in 
a period range of 0.4 - 10 day s. The upper limit was chosen 
because lAigrain et all [|2007h show that the time sampling 
of our observations, combined with the geometric probabil- 
ity of a transit, is insensitive to planets with orbital periods 
greater than ~10 days. The lower limit was selected below 
the typical boundary of ~1 day because we wanted to search 
for systems with extremely short periods, and because we 
wanted to test our sensitivity at these short periods in the 
simulations described in Section 5. It is important to remem- 
ber, however, that there are no confi rmed planet dete ctions 
to date with periods less than 1 day. (|Sahu et al.ll200fj found 
5 planetary candidates with periods less than a day, and as 
low as 0.42 days, however, they remain unconfirmed because 
they are too faint for follow-up spectroscopy.) We show our 
transiting system recovery fraction as a function of period 
and the total number of required transits for a positive de- 
tection in Figure [7] The figure gives the probability of 1, 2, 
or 3 transits being present in the data as a function of orbital 
period. Transits are considered to be present if there are data 
within the range of phases <j> < 0.1W/P or <f> > 1 — 0.1W/P 
where W is the expected transit width and P is the orbital 
period, i.e. a transit is present if we observe some portion of 
the central 20% of the transit. Given that ^ 3 detections are 
needed to accurately the determine the period of a system, 
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Figure 7. Fraction of transiting systems recovered as a function 
of orbital period. The results for 1 (top, solid curve), 2 (middle, 
dotted curve), and 3 (bottom, dashed curve) detected transits are 
shown. 



Figure [Jj shows that we can only hope to recover reliable 
periods for systems with orbital periods < 2 days. 

While searching each light curve for the best-fitting oc- 
cultation model, we also fit for the best brightening model 
(defined as an increase in brightness, rather than decrease 
as is the case for an actual occultation) as described in 
iBurke et al.l l|2006l ). This secondary fit has little impact on 
the numerical efficiency of the fit, because the algorithm si- 
multaneously searches for the greatest reduction in x 2 in 
terms of both brightening and dimming. We have no reason 
to believe that objects with correlated noise on the time scale 
of an occultation or non-periodic variable objects not se- 
lected by the method described in Section 4.1 preferentially 
show correlated decreases, rather than increases, in flux. On 
the contrary, one would expect that for these cases both 
an occultation and brightening could provide good models 
to the data. Therefore, we use the ratio of the improve- 
ment in x 2 f° r an occultation model to a brightening model 
to discriminate against non-eclipsing objects by requiring 
Ax 2 /Ax?- > 3, where Ax 2 is the improvement in x 2 using 
an occultation model and Ax- is the improvement in x 2 
relative to a brightening model. 

Figure|5]shows a plot of Ax 2 - against Ax 2 for every light 
curve in our sample. The solid line shows the boundary of 
our requirement Ax 2 /Ax~ > 3. Objects below the line pass 
the selection criteria and are modelled significantly better 
by an occultation than by brigh tening. 

We employed the method of I Pont et alj l|2006l ) to deter- 
mine the significance of a detection for each light curve in the 
presence of red noise. They show that the detection statis- 
tic, Srcd in a light curve with red noise can be found without 
a fit to the individual white and red noise components. We 
briefly summarise the procedure here: first we found the best 
fit solution assuming white noise (this is equivalent to max- 
imising S), then we masked points in-occultation. We then 
calculated the mean flux over a sliding interval equal to the 
duration of the detected occultation, where the sliding steps 
are smaller than the interval between flux measurements. 
We grouped those flux measurements into bins based on the 
number of data points in the sliding interval, and we cal- 
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Figure 8. The improvement in x 2 with respect to a brighten- 
ing model (Ax 2 ,) as a function of the improvement in \ 2 with 
respect to an occultation model (Ax 2 ) for every searched light 
curve. The solid line corresponds to A% 2 /Ax^ = 3. We require 
points to lie below the line to be considered occultation candi- 
dates. The green diamonds indicate light curves that pass all of 
our selection criteria. The blue X's are a subset of the green dia- 
monds. In addition to passing all our selection criteria, the blue 
X's have observed occultation depths that may be compatible 
with a planetary companion (see Section4.4). 



15 16 17 18 19 

i-mag (detector) 

Figure 9. Detection statistic accounting for red noise as a func- 
tion of CCD i-band magnitude for all candidate cluster members 
brighter than I = 19 (black dots). Objects with large values of 
•S'rcdi which occur over the full magnitude range of our observa- 
tions, indicate our best occultation candidates. The red triangles 
show stars classified as rotators. The dashed horizontal line shows 
our adopted threshold in S Ie d- Stars above the line meet the cri- 
teria of the threshold. The green diamonds and blue X's are the 
same as Figure [8] 



culated the variance of the flux measure ments in ea c h bin . 
These variances give an estimate to what I Pont et al.l (|2006t ) 
call the V(n) function, where V(l) equals the variance in 
bins with only one point, V(2) equals the variance in bins 
with two points, and so on. The detecti on statistic can b e 
measured with V(n) using Eqn. 7 from I Pont et al.l (2006), 
reproduced here: 



d A 



(4) 



where d is the depth, n is the total number of data points 
in-occultation, the sum is over all occultations, and nt is the 
number of points in the fcth occultation. 

After measuring S 2 for each light curve, we determined 
the corresponding value of S TB d- We show S Te d as a function 
of detector i-band magnitude in Figure [9] for every object 
in our sample. From the figure it can be seen that we found 
strong detections (S Ic d > 6) over the full magnitude range 
of our observations. The fact that there is no bulk trend 
in Srod as a function of magnitude shows that S Ic d takes 
into account any residual correlated noise, which is typically 
magnitude dependent. The triangles show objects classified 
as rotators following the procedure in Section 4.1. Figure 
indicates that our filtering method is successful, as the vast 
majority of the rotators have detection statistics consistent 
with non-eclipsing light curves. A plot of S Ic d as a function 
of best fit period can be seen in Figure 1101 We find a large 
over-density of points at a period of 0.5 days. A histogram of 
the best fit period for every object shows a high frequency at 
0.5 and 1 days. Visual inspection of these light curves shows 
a systematic effect where slight changes in the median flux 
(~ 0.005 mag) from one observing season to another can 
be fit extremely well with periods of ~0.5 or ~1 days. We 
are not entirely certain why the flux changes from season to 
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Figure 10. function of best fit occultation period for 

all candidate cluster members in our sample. The vertical lines 
show regions of detected periods that are excluded to remove 
false positives, as described in the text. The dashed horizontal 
line shows our adopted threshold in S TC d- Stars above the line 
meet the criteria of the threshold. We note that the scale of S re d 
has been reduced, which excludes the two candidate occultation 
systems at (i, 5 rcd ) « (17.1, 40) and (17.6,25) in Figure|9] The 
green diamonds and blue X's are the same as Figure [8] 



season for this subset of objects, but we wish to eliminate 
these false positives. Therefore, we do not consider any light 
curves with a best fit period within 1 ± 0.05 and 0.5 ± 0.01 
days as detections. The vertical lines in Figure \W\ indicate 
these regions of discarded periods. 
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4-3.1 Detection Threshold 

Following their examination of the effects of corr elated noise 
on the ability of transit surveys to detect planets, IPont et al.l 
(2006) use simulated light curves to show that false positives 
rarely have values of 5 re d larger than 7. They then argue 
that the adopted detection threshold will vary from survey 
to survey, but that it will typically be in the range 7-9. 

Following from this, we chose to adopt an inclusive de- 
tection threshold of Sdet = 6.5. Figures [9l and [101 show that 
very few light curves have S ro d > Sdet, while a histogram 
of S'red for all objects within the allowed range of best fit 
occultation periods shows very few counts when S ro d > 6.5. 
We acknowledge that the selection of this cut will likely in- 
crease our false positive rate relative to other transit surveys, 
however, in our case this may not prove to be entirely detri- 
mental. We investigate the choice of our detection threshold 
in Section 5, and find that our Monte Carlo simulations cor- 
roborate our adopted detection thresholds. 

There are a few circumstances that make our situation 
fairly unique. The first is t hat the size of the clu ster matches 
the FLAMES field of view (|Aigrain et al.ll2007T l. The second 
is the high scientific value of any transit candidates which 
turn out to be cluster EBs. Therefore, we can afford to adopt 
a relatively low threshold and accept more false positives 
than a typical survey, because (1) we can follow-up many 
of these candidates at the same time and (2) the detection 
of any occulting body in the cluster provides an important 
discovery. In Section 5, we use a series of Monte Carlo sim- 
ulations to estimate our sensitivity and false positive rate 
based on our selection of Sdet- 



4.4 Occultation Candidates 

15 out of the 1180 candidate cluster members pass all of our 
selection criteria as occultation candidates. Following depth 
considerations, planetary transits can be ruled out for all but 
six of the candidates. For the cases where transits cannot be 
ruled out very large planets, > 1.5i?j up , would be needed to 
explain the observed occultation depth. 

We estimate the occultation depth for each of the 15 
candidates in order to exclude any systems that could not 
be a dwarf-planet system from our list of transit candidates. 
The depth is measured by eye in the following way: the min- 
imum flux during occultation is subtracted from the out of 
occultation flux on the same night. This results in an esti- 
mate of the depth in magnitudes. We measure the occulta- 
tion depth by eye because intrinsic stellar variability makes 
the best fit depth from our occultation search unreliable. 
We do not rule out planets in any candidate where the ob- 
served occultation depth could be explained by a planet with 
R p ^ 3-Rjup. This leaves six candidates as possible transits, 
while the other nine systems remain EB candidates. We de- 
fer a discussion of the EB candidates to a later paper. In 
Table [2] we include the observed occultation depths and cor- 
responding minimum planet radii necessary to explain the 
occultation. We note that the uncertainty in our estimates 
of the necessary planetary radius, which relies on stellar evo- 
lution models and a visual measurement of the occultation 
depth, is large. Table[5]also contains the occultation param- 
eters of each of our transit candidates. AxV^X- and S'red 
are the detection statistics discussed in Section 4.3. 
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Table 2. Fit parameters for candidate transiting systems in NGC 
2362. Star is the object identification number from this work. S 
is the occultation depth, measured by eye, and R p is the corre- 
sponding planet radius necessary to account for a transit of that 
depth assuming a central transit. 

Table[3]summarises the properties of the six stars where 
we cannot rule out the possibility of a transit. We include 
the survey identification number, J2000 coordinates, optical 
photometry taken at CTIO, NIR photometry from 2MASS 
where available, as well as the model mass and radius for 
each candidate from our adopted magnitude- mass-radius re- 
lation (see Section 3) . These estimates of the stellar parame- 
ters assume that the observed flux is from a single star, and 
that there is no occulting body which contributes to the to- 
tal brightness of the system. If these systems are EBs, then 
these values serve as upper limits to the mass and radius of 
the primary star. 

In Figure[TT]we show the full light curves for each of our 
transit candidates, as well as a simulated transit included as 
a reference. We elect to show the entire light curve because 
in most cases, with the exception of a few, we do not detect 
enough occultations to accurately determine the period of 
the system. We show the observed light curves before they 
have been filtered for variability, because out of occultation 
variability is indicative of youth, and hence cluster member- 
ship. Note that for the model transit the duration is shorter, 
and the depth is shallower than for our candidate systems. 
This indicates that there is a good chance that our transit 
candidates are in fact EB systems. 

We note that for each case multi-epoch mid- to high- 
resolution spectroscopy is going to be needed for RV mea- 
surements in order to (1) determine cluster membership via 
the systemic velocity and (2) confirm whether the observed 
reductions in flux are the result of an occulting planet or 
dwarf. Toward these two goals, we will observe NGC 2362 
with VLT/FLAMES. In most of these cases additional pho- 
tometric observations will be needed in order to provide bet- 
ter orbital phase coverage as well as additional detections of 
the single or partial occultations already observed. 

5 RESULTS 
5.1 Simulations 

We report here on the Monte Carlo simulations designed 
to characterise the sensitivity of our survey to transits and 
estimate the false positive rate. We do this by inserting tran- 
sits into the raw data, correcting for correlations with seeing 
(see Section 3), and searching the resulting light curves us- 
ing the same procedure described in Section 4. We consider 
all simulated light curves to be positive detections if prior to 
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Table 3. Candidate transiting systems in NGC 2362. Star is the object identification number from this work, a and S are the RA and 
DEC, respecitively. V, V — R, and V — I are optical photometry measurements from this work. J, J — H and H — K are IR measurement 
from 2MASS, where available. M* and R* are the model mass and radius, respectively, using our adopted magnitude-mass-radius relation, 
and assuming the candidate is being occulted an object or material which does not contribute to the total brightness of the system. 
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Figure 11. Full light curves for the six stars where we could 
not rule out the planetary hypothesis. From top to bottom: 3- 
10048, 3-3739, 3-7559, 5-6469, 6-9484, 7-4723. The bottom panel 
shows a light curve with an inserted transit (see Section 5). The 
transits have been highlighted relative to the out-of-transit data. 
Note that the simulated transit is typically shorter in duration 
and shallower than the events shown above. Gaps in our observa- 
tions have been removed and are indicated by the vertical lines. 
For reference the relative night of observation is listed above the 
corresponding portion of the light curve. 



the injection of a transit A^-/Ax+ < 3 and/or S TC d < 6.5, 
whereas following the insertion of a planet Ax 2 - / Ax+ > 3 
and Srod > 6.5. Saturated stars and stars fainter than / = 19 
were not included in the simulations. We only insert transits 
into 1171 of the 1180 light curves we initially searched for 
occultations in Section 4. The 9 stars where we ruled out 
the planetary hypothesis based on the depth of the occulta- 
tions were excluded from the simulations, because they have 
large reductions in flux that would result in positive detec- 
tions regardless of the shape or size of the inserted transit 
signal. 



We use the simulations to characterise the survey's sen- 
sitivity to HJs. We define the sensitivity, 5, as 



5 = 



rid 

N, sirr. 



(•>) 



where rid is the total number of planetary systems that we 
detect and iV s ; m is the total number of simulated planetary 
systems. This represents a measure of our ability to detect 
planets over the entire range of possible inclinations, not 
just those planets that transit their host star. From S, we 
can determine the expected number of planet detections, 
following some assumptions about the fraction of stars with 
planets. 

These simulations accurately address the actual signal 
to noise in our data, the real cadence and observing win- 
dows in our data, and, perhaps most importantly, the intrin- 
sic variability present in the light curve of each individual 
star. This rep r esent s an improvement to the simulations in 
lAigrain et all (|2007h . which assume a signal to noise based 
on the average noise properties of non-variable light curves 
and ignore the intrinsic variability of individual stars. 

When testing the sensitivity of a transit su rvey it is im- 
porta nt to inject realistic transits into the d ata jBurke et aL 
2006). Therefore, we adopt the formalism of Mandcl fc Ago] 
(2002), which provides a method for calculating analytic 
transits given the ratio of planet and star radii as well 
as the geometry of the orbit. Specifically, we use the 
quadratic limb- d arken ing model presented in Section 4 of 
iMandel fc Ago! (|2002l ). We determine the limb-darkening 
coefficients from the tables of IClaretl (|2000T ), and the stel- 
lar properties come from the isochrones of the NextGen, 
DUSTY, and COND models, as described in Section 3. We 
acknowledge the perils of assuming the validity of the very 
models we are trying to test, however, without the use of 
some model to relate absolute magnitude to stellar mass 
and radius the following simulations would be impossible. 

For our simulations we assume circular orbits. Given 
the relatively low eccentricities of all discovered extrasolar 
planets orbiting within 0.1 AU, this is a reasonable assump- 
tion. For planets with semi-major axis less than 0.1 AU the 
median eccentricity is 0.04, with a mean eccentricity < e > 
= 0.08, and only two planets with e > 0.30 

Each simulation was run using a uniform distribution 
in logio period from 0.4 to 10 days, while the epoch was 
selected randomly from a uniform distribution of orbital 

4 http://exoplanet.eu/catalog.php 
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phase. Finally we selected an orbital inclination from a 
uniform distribution in sini. Transits only occur when the 
planet passes in front of the stellar disc, i.e. only in cases 
where 



(.R* +R p )/a ^ cosi, 



(6) 



where R* is the stellar radius, R p is the planet radius, and a 
and i are the orbital semi-major axis and inclination, respec- 
tively. Therefore, we could choose to select to only simulate 
the cases where the inserted planet transits the stellar disc. 
We choose to select i from the full range of orbital incli- 
nations, however, because we want to use the simulations 
to test our sensitivity to all planets, regardless of whether 
the planet transits its host star or not. This choice has a 
negligible effect on the cost of the total computing time. In 
the cases where the planet does not transit (i.e. i does not 
satisfy Eqn. |6}, or the phase coverage is such that we do 
not observe any transits, we do not need to run the search 
procedure as the light curve remains unchanged and the re- 
sults of the occultation search are the same as those already 
obtained during the initial search procedure described in 
section Section 4. 

We inserted planets of radii li?j up and 1.5-Rj up into the 
light curves. We arrived at the lower planetary radius be- 
cause ~ 1-Rjup is roughly the lower limit of rad ii we are sensi- 
tive to following the preliminary simulations of lAigrain et al.l 
120071 ) . W e chose the large r radiu s based on the evolutionary 
models of lBurrows~ et al. ( 1993) for ex trasolar giant plan- 
ets. According to Burrows et al.l ( |l997i ) at an age of 5 Myr 



planets in the mass range 0.7Mj up to 2Mj up have radius 
R p ~ 1.5-Rj up , while the largest radius occurs for a Saturn 
mass planet (0.3Mj up ) with R p ~ 1.7iijup.0 

For every star included in the simulations we indepen- 
dently inserted 300 planets with R p — li?j up and 300 planets 
with R p — 1.5ii.ju P . Specifically, for each iteration of the sim- 
ulations we would: (1) add a planet to the system with given 
R p , i, and a, which comes from the adopted stellar mass and 
randomly selected period, (2) check to see if Eqn. \§\ is sat- 
isfied, if it is then (3) insert the transit into the light curve, 
and (4) check to see if the randomly chosen epoch results 
in any observed transits. In total this corresponds to the in- 
sertion of 708 000 simulated planetary systems of which we 
searched ~ 101 500 and ~105 500 systems with observed 
transits for the inserted li?j up and 1.57?j up planets, respec- 
tively. 



5.1.1 False Positives 

False positives cannot be identified by a light curve alone. 
RV measurements, which can determine the planetary na- 
ture and cluster membership of any candidates, are needed 
to fully characterise the fra c tion o f light curves that are in 
fact false positives. iBrownl (2003) calculates the expected 
number of false positives due to stellar companions for tran- 
sit survey s of the field s tar population. An analysis similar 
to that of IBrownl |2003h adapted to the stellar population 
in NGC 2362, which is different from field stars, would be 
desirable, but is considered beyond the scope of this paper. 
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2.0 
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0.025 
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0.025 


5.0 


2.0 


0.066 


0.022 


0.132 


0.022 


5.5 


2.0 


0.062 


0.020 


0.128 


0.020 


6.0 


2.5 


0.048 


0.011 


0.117 


0.011 


6.5 


3.0 


0.038 


0.005 


0.107 


0.005 


8.0 


3.0 


0.029 


0.003 


0.095 


0.003 


10.0 


5.0 


0.012 


0.002 


0.072 


0.002 


15.0 


10.0 


0.001 


0.000 


0.038 


0.000 



Table 4. Sensitivity, S, assuming a logarithmic distribution in 
orbital periods, and the upper bound to the fraction of false posi- 
tives {TV) as a function of planet radius. Note that S and TV are 
not directly comparable because the quoted values of S assume 
every star has a HJ. 



5.1.2 Sensitivity as a Function of Detection Threshold 

We now evaluate our selected thresholds based on the results 
of our simulations. In Table [4] we summarise the sensitivity 
and maximum false positive fraction for a number of differ- 
ent thresholds in S rB d and Ax 2 /Ax 2 _ for both li?.j up and 
1.5i?j up . As noted above, we cannot determine the number 
of false positives from the light curves alone. If we assume 
that there are no transiting planets in our data, then we can, 
however, determine an upper bound to the number of false 
positives. This upper bound is determined by assuming that 
all systems which pass the detection criteria are considered 
false positives: 



TV = 



^cand 

1171 ' 



(7) 



http : //zenith . as . arizona . edu/>~ burrows/ 



where TV is the upper bound on the number of false pos- 
itives, n can( j is the number of systems that pass all the se- 
lection criteria, and 1171 is the total number of searched 
systems. 

The quoted values of S assume that every star has a 
HJ, which means the number of detections relative to the 
number of false positives cannot be found by dividing S by 
TV . To find this ratio S would have t o be multiplied by the 
fraction of stars with a HJ, or ~1% l|Gaudi et al J 1200a ; we 
caution that this result is limited to the solar neighbourhood 
and may not apply to a cluster at the distance of NGC 
2362, where the formation environment may be significantly 
different than that in the solar neighbourhood). 

What Table [4] clearly shows is that there is very little 
improvement in the sensitivity when the S te d threshold is 
lowered below 5.5, while at the same time TV grows rel- 
atively quickly. At the same time, to obtain TV < 0.1% 
would require the selection of very large detection thresh- 
olds. Given the relatively small number of cluster stars in 
this survey, such a large threshold would be inadvisable. 
Reducing the thresholds from these large values to S rc a ~ 
6 and Ax 2 /Ax- ~ 3 results in the largest gains in <S rela- 
tive to the increase in TV. We arrive at our final thresholds 
of Sred = 6.5 and Ax 2 /Ax 2 _ = 3, because reducing each of 
those thresholds by 0.5 creates additional occultation can- 
didates whose variability is clearly not the result of eclipses 
or transits. 

If we assume that 1% of stars have a HJ then the values 
quoted in Table [4] are somewhat discouraging. This would 
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mass 
(M ) 


radius 

(%) 


l.Ofljup 


1-5-Rjup 


S 


5 


< 0.3 


< 0.92 


0.048 


0.114 


0.3 - 0.5 


0.92 - 1.18 


0.038 


0.104 


0.5 - 0.7 


1.18 - 1.38 


0.027 


0.102 


> 0.7 


> 1.38 


0.037 


0.109 



Table 5. Sensitivity as a function of stellar mass and radius. The 
mass and radius are determined from our adopted magnitude- 
mass-radius relation as described in the text. S is the sensitivity 
for each range of masses. 

mean that, given our adopted thresholds, for a large sur- 
vey in order to detect 1 planet at 1/Zj up and 1.5J?j up wc 
would have ~13 and ~5 false positives, respectively. This 
is somewhat abated considering this cluster lies within a 
singe FLAMES field of view, meaning that the telescope 
time needed for spectroscopic follow-up is essentially inde- 
pendent of the number of candidates. 



5.1.3 Sensitivity as a Function of Stellar Mass 

As noted previously, the presence of a planet cannot be in- 
ferred from a light curve alone. Given that consistent RV 
variations are needed to confirm a planet, it would be very 
useful to understand our sensitivity as a function of stellar 
mass for stars in the cluster. Our ability to detect planets 
in NGC 2362 via RV vari ations is strongly de pendent on 
the mass of the host star (|Aigrain et al.|]2007t ) . Assuming 
a Jupiter mass planet, the greatest signal in RV variations 
is going to occur for the lowest mass stars, however, these 
stars are also going to be the faintest in the cluster meaning 
it will be difficult to obtain spectra with sufficient signal-to- 
noise to detect the RV variations. While moving to higher 
mass stars will dramatically improve the signal-to-noise in 
a single spectrum it will cause a significant reduction in the 
RV amplitude, again assuming a Jupiter mass companion. 
This will make it difficult to detect a planet. Thus, an un- 
derstanding of the sensitivity as a function of stellar mass 
becomes extremely important when determining which sys- 
tems to follow-up for this cluster and for the Monitor project 
as a whole. 

In Table[5]we summarise the sensitivity to HJs as a func- 
tion of stellar mass. We also show the corresponding radii 
from our adopted magnitude-mass-radius relation. Some- 
what surprisingly, Table [5] shows that we are more likely 
to recover a planet around a smaller star. This is surprising 
because the stars with the lowest pre-transit rms are all at 
the bright (and hence higher mass) end of the cluster. This 
implies that the effect of reduced stellar radii is more impor- 
tant than a low rms when t rying to detect transits. This con- 
firms the initial findings of lAigrain fc Pontl (|2007l ). who ex- 
amined the dete c tabili ty of transits in a hypothetical cluster. 
lAigrain fc Pont) i|2007l ) found that when red noise is consid- 
ered it becomes more difficult to detect planets around the 
brighter cluster members even though they have a smaller 
rms. Our results suggest that small stellar radii may be the 
most important factor in detecting transits given that our 
survey is most sensitive to planets orbiting stars with mass 
< O.3M and radius < O.92R . 



5.2 Expected Number of Hot Jupiter Detections 

Using the Monte Carlo simulations described above we esti- 
mate the expected number of detectable short period planets 
in NGC 2362. Many factors must be accounted for in order 
to estimate the number of expected detections, including: 
the frequency of HJs, the total number of candidate cluster 
members and the contamination from field stars, the geo- 
metric transit probability as a function of period, and our 
recovery rate as a function of period. We note that the rela- 
tively small number of cluster members in NGC 2362 means 
that a null result in NGC 2362 does not carry much signifi- 
cance. 

We begin by assuming that the frequency of large, short 
period pla nets is the same as that found in the solar neigh- 
bourhood. iGaudi et al.l l|2005l ) identify two empirically de- 
fined populations of planets with periods < 10 days. Follow- 
ing from their work we assume the frequency of hot Jupiters, 
O p i, with periods of 3-10 days, t o be 1 %. Strictly speaking, 
the estimates from lGaudi et all (j2005T) include only planets 
with periods up to 9 days. We extend their estimate to 10 
days without loss of generality. The frequency of HJs in the 
other population with periods < 3 days, which [Gaudi et all 
l|2005l) refer to as VHJs, is ~0.15%. This is then multiplied 
by the number of candidate cluster members, JV», and the 
probability of membership, P m emb, to arrive at the expected 
number of HJs in our data set. The number of HJs is not ran- 
domly distributed in orbital period. The actual distribution 
of orbital periods is currently ill constrained, though it likely 
depends on the physical parameters of the system in which 
the planet is formed. In order to calculate the number of 
expected detections we must make some assumption about 
the period distribution. Therefore, we assume a uniform log- 
arithmic distribution in orbital period. We then separate the 
expected number of HJs into equal size bins in log(period) 
space and calculate the number of planets in each period bin 
i, which we then sum to arrive at 

n 

iVdet = 5^iV P i(perO, (8) 

i — m 

with 

tO pl f pk S t (9) 

where iVdct is the total number of detectable HJs present in 
the data set, N*P mem \, is the number of cluster members we 
observed, O p i is the frequency of HJs in the solar neighbour- 
hood, / p ij is an estimate of the fraction of the total number 
of HJs in bin i based on a logarithmic distribution of orbital 
periods, and Si is our sensitivity to planets in the ith pe- 
riod bin. Our final assumption is that there are no HJs with 
periods less than 1 day, because to date there have been no 
RV confirmed planets found with periods < 1 day. 

The resulting expected distribution of HJs is shown in 
Table [S] We find that we would only expect to detect 0.19 
HJs in NGC 2362, which is consistent with a null detec- 
tion. The values quoted in Table [6] assume that HJs orbiting 
young stars have bloated radii of 1.5J?j up . If instead we as- 
sume that these planets all have radii of 1.0i?j up then the 
sensitivity is reduced and the expected number of detectable 
HJs in our data is reduced to 0.038. Under either circum- 
stance, it is clear that we would not expect to find any HJs 
in our data. 
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1.5 


6.61 


3.0 - 10.0 




0.453 


0.697 



Table 6. Number of expected detectable, transiting planets in 
NGC 2362. We have separated the results into 3 different popula- 
tions based on orbital period. Fp\. is the expected number of HJs 
in our data in the corresponding period range. Si is the sensitivity 
to 1.5-Rjup HJs, as determined by our Monte Carlo simulations. 
Npi(Per) is the expected number of transit detections at each 
period. 



5.2.1 Upper Limits on the Planetary Fraction in NGC 
2362 

Given our null detection (assuming that none of our candi- 
dates are planetary transits, which has yet to be confirmed) 
we can place upper limits on the fraction of stars in NGC 
2362 with a HJ. We do this by assuming that all stars in 
the cluster have a HJ, and then calculate the correspond- 
ing number of expected detections based on the results of 
our Monte Carlo simulations. Following from there we can 
evaluate the statistical significance of a null result given the 
number of expected detections. 

Based on RV surveys of the s olar neighbourhoo d, we 
know that not every star has a HJ (jGaudi et alj|200a ). The 
actual number of stars with detectable HJs will be a Poisson 
distribution with mean n = fpNdet, where f p is the actual 
fraction of stars with HJs and N de t is the expected number 
of detections assuming all stars have a HJ. For a Poisson 
distribution the probability of n detections given a mean fi 
is 



(10) 



For our null result n = 0, which when we substitute into 
Equation [10] yields 



P(0;f p N dBt ) = e ~ / " Jv ^. 



(11) 



To obtain an upper limit on f p at significance a we require 
that 



a P(0;f p N det ) 



(12) 



When we substitute Equation [TT] into Equation[l2]we arrive 
at 

In a 



(13) 



which allows us to place an upper limit on f p at any sig- 
nificance q (or confidence level 1-a). For example, if one 
expects to detect 3 planets, then there is a 5% chance of de- 
tecting zero planets from Poisson statistics alone (Eqn. [TT|) . 
Equivalently, to convert a null detection to a statistical re- 
sult at 95% confidence or greater, the expectation value of 
the number of planets detected should be 3 or greater. 

In Table [7] we show the derived upper limits on f p for 
q = 0.05 and 0.01 (corresponding to confidence levels of 
95% and 99%, respectively) for HJs with R p — 1.0i?j up 
and 1.5-Rjup. We quote upper limits for two different period 
ranges: 1-3 d, and 3-10 d. In order to calculate the upper 



Table 7. Upper limits on the fraction of stars with HJs. R p is 
the planet radius. N dct is the expected number of detections if 
every star has a HJ. f p is the fraction of stars with HJs in the 
given period bin, and a is the significance level of the upper limit 
on fp. *For the case of 3 - 10 day 1-Rj up planets we cannot place 
upper limits below every star in the cluster having a HJ in this 
period range. 



limits in each period bin we assume that every star has one 
planet of radius R p . We then determine the number of stars 
with a HJ in a given bin, assuming a uniform logarithmic 
distribution, and multiply this by the sensitivity, Si, in the 
given bin to arrive at A/det- The most significant upper limits 
we find are at a = 0.01, where if HJs at ~5 Myr are about 
1.5-Rjup then the upper limit on short (1-3 days) period HJs 
is 22% while the upper limit on HJs with periods between 3- 
10 days is 70%. These limits provide strong evidence against 
the 'survival of the lucky few' scenario, where most planets 
migrate into their host star and only those that form just 
prior to the d ispers al of the disc survive, as discussed in 
lAigrain et al.l (|2007l ). We note that if HJs at this age tend 
to be closer to 17?j U p we cannot place reliable constraints on 

5.3 Discussion 

Given the number of cluster stars (~475), the number of 
hours observing (~100), and the faintness of many of the 
cluster members it is unlikely that we would have detected 
a planet in NGC 2362. The failure to detect a planet is un- 
surprising based on our simulations, and the simulations of 
lAigrain et al.l ([2007) , which predicted zero detectable plan- 
ets in our data for this cluster (again, we are assuming that 
none of our transit candidates are in fact transiting planets, 
which remains to be confirmed). In fact, the detection of a 
planet would have been more inconsistent with the expecta- 
tions than a non-detection, and would likely provide strong 
evidence for a greater incidence of planets around young 
stars than MS stars. 



5.3.1 Comparison with Other Cluster Surveys 

There have been many surveys searchin g for planetary tran- 
sits in star clusters (see Wcldrakc 2007 for a review). These 
surveys are unique relative to the many shallow, wide-field 
transit surveys, because they examine stars of a known age 
and metallicity. Therefore any planets or low mass EBs 
found in these clusters can be used as observational con- 
straints on stellar environments and their evolution. 

These cluster surveys also present an opportunity to 
place upper limits on HJ incidence at a number of stellar 
ages. Unfortunately, few of these surveys have actually cal- 
culated upper limits o n fp, but we can co mpare our results 
with those that have. iBurke et al.l {2006) surveyed the ~1 
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Gyr cluster NGC 1245, and found 95% confidence upper 
limits for 1-3 day orbits which are a factor of ~2 smaller 
than the upper limits from this work. They find upper limits 
of 6.4% and 24% for 1.5i?j up and 1.0J?j up HJs, respectively, 
compared to upper limits of 14% and 52%, respectively, from 
this work. We report an upper li mit for 3-10 day peri od HJs 
of 45%, which is lower than the [Burke et all i|2006l ) value 
of 52% for 1.5i?.jup planets. Neither study was able to place 
meaningful upper li mits on the incide n ce of 1-Rj up planets 
in 3-10 day orbits. [Bramich fc Hornel (120061 ) were able to 
place better upper limits on f v than those found in this 
study, however, they included both field and c luster stars 
(N, ~ 30000) in their analysis. IWeldrake etafl (|2008bl) ob- 
served 31 000 stars in the globular cluster lu Centauri for 25 
nights. These observations allowed them place a 95% con- 
fidence level upper limit of 0.1% on 1.5J?j up planets in 1-3 
day orbits. Given the significan t ly lar ger nu mber of targets 
in both the iBramich fc Hornel (2006f ) and IWeldrake et all 
(2008b]) surveys it is unsurprising that they place more re- 
stricting upper limits on f p . 



5.3.2 Extensions to the Remainder of Monitor 

We can also explore what happens when we extrapolate our 
results to Monitor as a whole. Monitor will observe a total of 
nearly 15 000 young (< 200 Myr) stars, most of which have 
not yet reached the MS. This significant increase over the 
number of targets in NGC 2362 will lead to large reductions 
in the limits on f p for young stars. 

In fact, observations of an additional Mo nitor target, h 
&y Per, which has ~7500 cluster members (|Aigrain et al.l 
2007), will lead to a considerable reduction in the limits on 
f p for young stars. Assuming the only difference between 
clusters is the number of observed star^f], then we would be 
able to reduce all the upper limits in Table [7] by a factor 
of ~13 following the observation of h & x P er - When we 
consider all the stars to be observed by Monitor we will be 
able to reduce the limits on f p by a factor of ~26. For the 
case of 1.5-Rj up planets in 3-10 day orbits this would mean 
an upper limit of ~2.7%. 

The above predictions assume that our achieved sensi- 
tivity in NGC 2362 will be the same in each of the clusters 
we have observed. Most of the Monitor targets are both older 
and less distant than NGC 2362, meaning in these cases we 
are probing stars with smaller radii and less intrinsic vari- 
ability. As seen in Section 5.1.3, smaller radii should result 
in an improved sensitivity. Assuming a sensitivity equal to 
that in NGC 2362 we can estimate the final upper limits on 
f p we would expect from Monitor. These values are sum- 
marised in Table [5] 



6 We note that this assumption is a significant over simplifica- 
tion of the actual situation given that the clusters are at different 
distances, are being observed with different telescopes, have dif- 
ferent noise properties, and are different ages, however, we pro- 
ceed simply to provide an order of magnitude estimate for the 
improvement in f p . 



R p 


per 


f P 
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(S = 52362) 


1.0 


1-3 


0.031 


1.0 


3-10 


0.183 


1.5 


1-3 


0.008 


1.5 


3-10 


0.027 



Table 8. Summary of the expected a = 0.01 upper limits on the 
fraction of young stars with HJs following ~100 hrs of observa- 
tions of each of the Monitor targets. R p and per are the planet 
radius and orbital period, respectively. f p is the upper limit on 
the fraction of stars with HJs assuming that we achieve a sen- 
sitivity in all other clusters equal to that which we achieved in 
NGC 2362. 

5.3.3 Future Observational Considerations 

Finally, we would like to discuss the limitations of our sur- 
vey of NGC 2362. We have learned that our sensitivity is 
largely limited by intrinsic stellar variability and stars with 
large radii. Both of these limitations could be reduced by 
observing older clusters: there will be less intrinsic variabil- 
ity while the stellar radii will be smaller given that the stars 
have had more time to contract toward their MS radius. 

Another way to increase the sensitivity of our survey 
would be to decrease the point-to-point rms and the noise 
over time-scales equal to or longer than the duration of a 
transit. Given that the red and white noise properties of our 
data are dependent on a number of factors, including the de- 
tectors, sky noise, and observing conditions to name a few, 
there is no simple panacea for reducing the noise. However, 
there is one slight change in observing strategy which would 
generally reduce the noise (both red, which would have the 
greatest effect on our brightest targets, and white, which 
would have the greatest effect on our faintest targets) in the 
light curves: conducting all observations within a single ob- 
serving season. Even in many of the light curves with low 
X 2 with respect to a flat model we notice small fluctuations, 
of order 0.01 mag, in the median flux following large gaps 
in our observations. These shifts are real and not the result 
of systematics. Making all the observations in a single ob- 
serving season will remove this slight source of variability 
from our data and create a greater sensitivity to transits. 
We also note that this change in strategy would have the 
added benefit of significantly reducing the computing time 
necessary to adequately search the light curves for transits 
using our search algorithm. 



6 SUMMARY AND CONCLUSIONS 

We have conducted a search for occultations in NGC 2362. 
We observed the cluster on 18 nights from February 2005 to 
January 2006 with the Mosaic II imager on the 4m Blanco 
telescope at CTIO. We achieved an average cadence of ~6 
minutes, which is sufficient to perform differential photom- 
etry and search the data set for any transits. We used a 
V,V — I CMD to photometrically select 1813 candidate clus- 
ter members. 

Following the selection of candidate cluster members, 
we developed a systematic method for searching and iden- 
tifying occultations in our light curves. This method con- 
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sisted of two major steps: (1) the identification and removal 
of intrinsic stellar variability due to rotation, and (2) the 
search for transits usin g the occultation search algorithm of 
lAigrain fe Irwinl l|2004 ) modified to account for red noise. 
Following the removal of saturated stars and stars too faint 
for spectroscopic follow-up, we searched a total of 1180 stars 
for transits, of which about ~475 are expected to be cluster 
members according to our contamination estimate. 

Our search identified 15 light curves with reductions 
in flux that passed all of our detection criteria. Only six 
of these systems, however, have observed variability that 
would be compatible with a planetary companion based on 
the observed occultation depth. Some of these systems may 
be cluster EBs, which would help to provide important con- 
straints on the mass and radius, and by extension evolution, 
of PMS stars. 

Using a series of Monte Carlo simulations we predict 
the number of detectable HJs and find it is consistent with 
our null result. With 99% confidence we place a limit on 
the fraction of stars in NGC 2362 with 1-3 d period HJs 
at <22%, while we limit the fraction with 3-10 d period 
HJs at <70%, assuming a planetary radius of 1.5J?j up . We 
compare these limits with other cluster surveys and find that 
we are not as sensitive to transits as other surveys, because 
we have far fewer targets. The limits for NGC 2362 provide 
observational constraints on the fraction of stars with HJs 
at an age < 10 Myr. From the simulations we also know 
that our sensitivity to transits increases as the stellar radius 
decreases, a so mewhat non-intuitive result, which supports 
the findings of lAigrain fe Pond (|2007T l. 

Finally, we examine the prospects of the Monitor 
project as a whole. If we assume the same sensitivity is 
achieved for the entire Monitor survey, we will be able to 
place an upper limit on the number of young stars with a 
1.5-Rjup HJ at ~2.7%, assuming a null detection. 
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